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Within the minimal supersymmetric extension of the Standard Model, the mass of the light 
CP-even Higgs boson is computed to three-loop accuracy, taking into account the next-to-next-to- 
leading order effects from supersymmetric Quantum Chromodynamics. We consider two different 
scenarios for the mass hierarchies of the supersymmetric spectrum. Our numerical results amount 
to corrections of about 500 MeV which is of the same order as the experimental accuracy expected 
at the CERN Large Hadron Collider (LHC). 
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I. INTRODUCTION 

Supersymmetry is currently the most-studied exten- 
sion of the Standard Model (see, e.g., Ref. [H). It pro- 
vides solutions to some profound theoretical problems of 
the Standard Model: the fine tuning of the Higgs mass, 
the (non-)unification of gauge couplings, a mechanism 
for spontaneous symmetry breaking, and a Cold Dark 
Matter candidate. 

The minimal supersymmetric extension of the Stan- 
dard Model (MSSM) is based on a two-Higgs-doublet 
model (2HDM) with five physical Higgs bosons: two 
CP-even h/H, one CP-odd A (also named the "pseudo- 
scalar" Higgs), and two charged scalars H^. Each par- 
ticle of this 2HDM receives a SUSY partner of opposite 
spin-statistics, where left- and right-handed components 
of a Standard Model Dirac fermion are attributed with 
separate scalars /l/r which mix to the physical mass 
eigenstates f x i%- 

Compared to the Standard Model, the MSSM Higgs 
sector is described by two additional parameters, usu- 
ally chosen to be the pseudo-scalar mass Ma and the 
ratio of the vacuum expectation values of the two Higgs 
doublets, tan/3 = v%/v\. The masses of the other Higgs 
bosons are then fixed by SUSY constraints. In partic- 
ular, the mass of the light CP-even Higgs boson, Mh, 
is bounded from above. At tree- level, it is Mh < Mz- 
Radiative corrections to the Higgs pole masses raise this 
bound substantially to values that were inaccessible by 
LEP [3, [3) 0|. The large numerical impact is due to a 
contribution ~ cvtAf 2 ~ Af t 4 coming from top- and stop 



quark loops (M t is the top quark mass and y/dtt is pro- 
portional to the top Yukawa coupling) . 

The one-loop corrections to the Higgs pole masses are 
known without any approximations [a, [3 0> ■ They 
show that the bulk of the numerical effects can be ob- 
tained in the so-called effective-potential approach in the 
limit of vanishing external momentum. Motivated by this 
observation, all presumably relevant two-loop terms have 
since been evaluated in this approach (for reviews, see e.g. 
Refs. More 

recently there has been quite some ac- 



tivity in the context of the MSSM with complex parame- 
ters which can lead to sizeable effects (see, e.g., Ref. [11|). 
The two-loop results are implemented in the numerical 
programs FeynHiggs [12] and CPsuperH H3 . [l3 | using on- 
shell particle masses, and in SoftSusy |l5l |. SPheno [III], 
and Suspect [13] using DR parameters, that is, dimen- 
sional reduction with minimal subtraction. The influence 
of terms that go beyond the approximation of vanishing 
external momentum has been investigated in Ref. [Til ]. 

Based mostly on the renormalization scale and scheme 
dependence, the theoretical uncertainty on the prediction 
of the light Higgs boson mass Mh has been estimated to 
3-5 GeV [HI E3- This is to be compared with the ex- 
pected experimental uncertainty of a Higgs mass mea- 
surement at the LHC of the order of 100-200 MeV [1|. 
At an International Linear Collider, this goes even down 
to roughly 50 MeV [2l| . These numbers clearly show the 
need for three-loop corrections to the SUSY Higgs bosons 
masses in order to fully exploit the physics potential of 
these colliders. 

In fact, quite recently the leading and next-to-leading 
logarithmic terms in ln(MsusY/-^t) at three-loop level 
have been obtained, where MsuSY is the typical scale 
of SUSY particle masses [22}. In this letter, we want 
to present the first genuine three-loop calculation of the 
lightest Higgs boson mass, focusing on a few simplifying 
limiting cases for the sake of brevity. In particular, we 
consider effects of order at a 2 , keep only the leading terms 
~ Mf, and neglect all mixing effects in the stop sector. 
More general results and their detailed phenomenological 
impacts shall be deferred to a later publication. 



II. THE HIGGS BOSON MASS IN THE MSSM 

At tree-level, the mass matrix of the neutral, CP-even 
Higgs bosons h, H has the following form: 
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The diagonalization of M.% tTee gives the tree-level re- 
sult for Mh and Mg, and leads to the well-known bound 
Mh < Mz which is approached in the limit tan/3 — > oo. 

Quantum corrections to the Higgs boson masses are 
incorporated by evaluating the poles of the Higgs boson 
propagator at higher orders. As mentioned in the In- 
troduction, the numerically dominant contributions can 
be obtained in the approximation of zero external mo- 
mentum (see, e.g., Refs. 23]) which we will adopt in the 
following. Furthermore, we will only consider corrections 
of order cttct 2 .. Apart from the quark, squark, and gluino 
masses, there is another parameter with mass dimension, 
the trilinear coupling of the soft SUSY breaking terms, 
A t . Before renormalization, we express it through the 
stop masses Mj , M t ~ 2 , the stop mixing angle 9 t , and the 
bilinear Higgs parameter /^susy as follows: 

2M t A t = (Ml - Ml ) sin 26 t + 2A/ tA t S usY cot /J . (2) 

The mass matrix M.% is obtained from the quadratic 
terms in the Higgs boson potential constructed from the 
fields fa and fa. They are related to the physical Higgs 
mass eigenstates via a mixing angle a. Since fa does not 
couple directly to top quarks, it is convenient to perform 
the calculations of the Feynman diagrams in the (fa, fa) 
basis. 

Including higher order corrections, one obtains the 
Higgs boson mass matrix 



M 



H 



X If .tree 
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which again gives the physical Higgs boson masses upon 
diagonalization. The renormalized quantities , E^ 2 
and £^0;, are obtained from the self energies of the fields 
fa, fa, A, evaluated at zero external momentum, as well 
as from tadpole contributions of fa and fa (see, e.g., 
Ref. i). Let us remark that if one sets Mf = Mj 2 and 
A t = 0, and evaluates only the leading contribution ~ 
M t 4 , then only ^ and the matrix M 2 H — M 2 H trcc 
is diagonal. On the other hand, if we allow for non-zero 
A t , also and E^^ contribute in general. 

The calculation of is organized as follows: All 
Feynman diagrams are generated with QGRAF [2J|. In 
order to properly take into account the Majorana char- 
acter of the gluino, the output is subsequently manipu- 
lated by a PERL script which applies the rules given in 
Ref. [25|. The various diagram topologies are identified 
and transformed to FORM |26j with the help of q2e and 
exp 2^, HI- The program exp is also used in order to 
apply the asymptotic expansion (see, e.g., Ref. (29j) in 
the various mass hierarchies. The actual evaluation of 
the integrals is performed with the package MAT AD [30I ] . 
resulting in an expansion in d— 4 for each diagram, where 
d is the space-time dimension. The total number of three- 
loop diagrams amounts to about 16,000. 



At three-loop level we need to renormalize the top 
quark mass, the top squark mass, and the stop mixing 
angle at the two-loop order. In addition, the one-loop 
counterterm of the gluino mass is needed for the renor- 
malization of the two-loop expression. We implement 
Dimensional Reduction (DRED) with the help of the so- 
called e-scalars which appear for the first time at two 
loops. The renormalization of the e-scalar mass is per- 
formed in the on-shell scheme, requiring that the renor- 
malized mass is equal to zero. In the literature this is 
referred to as DR' scheme. 

The one-loop on-shell counterterms are well-known 
(see, e.g., Refs. [1, [U 32, H^)- As far as the two- loop 
counterterms for the squarks and quarks are concerned, 



one can find the results in Refs. [34j, [35(. However, it is 
rather tedious to extract the results for the mass hier- 
archies we are interested in. Thus, we re-computed the 
corresponding corrections. 

To our knowledge, the two-loop counterterm for the 
stop mixing angle is not yet available in the literature. It 
turns out that in our approximation, where Mj = Mf 2 
and At = 0, only the one- loop counterterm of the mixing 
angle enters the three-loop result. 

As a cross check for our calculation, we recalculated 
the exact two-loop result (in the limit of vanishing ex- 
ternal momentum) and find perfect agreement with the 
literature [23I [36l]. 

Furthermore, the expansion of the exact expressions 
confirms the limiting cases discussed below. Both the 
two- and three-loop calculations are performed for a gen- 
eral QCD gauge parameter £g. The independence of the 
final results on £5 serves as another welcome check on 
the correctness of our result. 

We use anti-commuting 75 which is allowed for fermion 
traces which involve an even number of 75 matrices. It 
turns out that all traces involving an odd number of 75 
vanish because they contain less than four gamma ma- 
trices. 

In the following we discuss three different cases for the 
mass hierarchy. In all cases we set the light quark masses 
to zero. 



(i) Super ■symmetric limit, i.e., Mt = Mj and the gluino 
and other squarks are massless: Mg — Mq = 0. The 
quantum corrections to the Higgs boson mass vanish in 
this case, as required by supersymmetry. Still, the in- 
dividual diagrams are different from zero and thus the 
calculation imposes a strong check on our setup. 

(ii) Massless gluino, M g = 0. Expanding in the limit 
M t <C Mi — Mq = Mstjsy, we obtain for the leading 
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with L Mt = ln(^ 2 /M 2 ) and L ts = ln(A/ 2 /Mf USY ). 

(iii) Common SUSY mass. In this scenario we assume 



M t < Mr 



Mr 



Msusy < Ms. Even though 



the top squark masses are equal and thus the mixing 
angle is zero, it is necessary to introduce a counterterm 
for 9t- Since the latter has contributions proportional to 
1/(M 2 — M? ), we expand the one- and two-loop result in 
this limit before inserting the counterterms. The cancel- 
lation of such terms in the final result provides another 
check on our calculation. 



It is important to keep A t ^ in the one- and two-loop 
contributions and to use Eq. ((2]) before renormalization, 
because the corresponding counterterms generate terms 
of order M t 4 at three-loop level. Setting A t = in the 
end, we obtain 
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where L t q = ln(M t 2 /M~). In Eq. © we only display the 
leading term in the 1/Mg expansion. We actually com- 
puted five expansion terms and observe a rapid conver- 
gence of the series — even for Msusy = Mj- It is inter- 
esting to mention that large cancellations occur among 
the cubic, quadratic, linear and non-logarithmic term of 
Yifa at three-loop order. E.g., for our default input values 
the sum of the cubic and quadratic logarithm is negative 
whereas the complete answer leads to a positive correc- 
tion for the a 2 coefficient of . 

If we express the result of Eq. ((5]) in terms of DR' 
parameters, we can compare with the results obtained 
22j. We find agreement both for the cubic and 



Rcf. 



quadratic logarithm. 



III. NUMERICAL RESULTS 



outlined in Ref. [38j which includes three-loop running 
and two-loop matching effects. As a result we ob- 
tain, e.g., a s (M t ) = 0.0926 for a common SUSY mass 
-^susy — 1 TcV. The SM input parameters are given as 
Gp = 1.16637 x 10" 5 GeV -2 , M z = 91.1876 GeV 
M t = 170.9 GeV 0. For the heavy squark mass (q / t) 
we use Ms = 2 TeV. 

In order to evaluate the tree-level approximation of 
the Higgs boson mass we also need the parameters M4 
and tan/3. If not stated otherwise we adopt the values 
Ma — 1 TeV and tan j3 = 40 . Since these parameters do 
not enter the corrections considered in this paper, they 
only have minor influence on the plots presented in the 
following. 

In Figs. [T]and [2] we discuss the difference between the 
Higgs boson mass evaluated with i-loop approximation 
and the tree- level result, 



In the remainder of this letter, we discuss the numer- 
ical effect of our result, restricting ourselves to At = 0. 
We adopt the on-shell scheme for the quark, squark and 
gluino masses. 

We choose fi — M t as the default value for the 
renormalization scale. First we compute a s (Mt) 1 de- 
fined in the DR scheme and the full SUSY theory, from 
the SM input value a s (M z ) = 0.1189 [13 which is 
given within five-flavour QCD. We follow the procedure 



AM 



(0 



M, 



(i — loop) 



M£ rcc . 
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Fig. [LJ shows AAff for i = 1 (dotted), i = 2 (dashed) 
and i = 3 (solid line) as a function of Msusy in the 
range between 200 GeV and 2 TeV. As is well known, 
the one- loop corrections are large, increasing Mh by up 
to 46 GeV. The two-loop effects are negative, reducing 
the size of the overall corrections by about 30% with 
respect to the one- loop result. 
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FIG. 1: AMh as a function of Msusy at one-, two-, and 
three- loop level. The renormalization scale is set to n — Mt- 



25.2 
25 

> 24.8 

o 

S 24.6 
< 24.4 
24.2 
24 



2- loop 

3- loop 

= 1000 GeV 



50 100 150 200 250 300 350 400 450 500 

u (GeV) 



IV. CONCLUSIONS 

To summarize, in this letter the three-loop corrections 
to the lightest Higgs boson mass have been computed 
in three different limits of the SUSY parameter space. 
For the phenomenologically interesting case where the 
gluino and top squarks have about the same mass and 
the remaining squarks are heavier, we observe effects of 
approximately 500 MeV. The dependence of the three- 
loop result on the renormalization scale indicates that 
the residual theoretical uncertainty matches the expected 
accuracy for a Higgs mass measurement at LHC and pos- 
sibly even at a future linear collider. 

It remains to say that the calculational setup which 
was used to obtain the results of this paper is not re- 
stricted to the specific MSSM parameter points consid- 
ered here. A more comprehensive study is in preparation 
and will be presented elsewhere. 
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FIG. 2: AMh as a function of the renormalization scale fj, at 
two- and three-loop level, where Msusy = 1 TeV has been 
chosen. 



The three-loop terms are much smaller and clearly sta- 
bilize the perturbative behaviour. At ji = M t , for exam- 
ple, they lead to a further reduction of AM^' by about 
400 MeV for Msusy = 300 GeV and an enhancement of 
about 500 MeV for M SUSY = 2 TeV. Note that the nu- 
merical impact is larger than the precision on the lightest 
Higgs boson mass as expected at the LHC. 

In order to estimate the size of the higher order cor- 
rections, we consider the dependence of the result on the 
choice of the renormalization scale. In Fig. [2] we plot 
AM^' as a function of fi which is varied from 50 GeV to 
500 GeV. The two-loop results show a variation of more 
than 1 GeV over this range. The error band derived in 
this way nicely covers the three-loop result, which itself 
varies by less than 35 MeV. For other values of Msusy 
the variation can reach up to 100 MeV. The three-loop 
curve in Fig. [2] shows a shallow minimum close to (i = M t 
which in turn is close to the intersection point of the two- 
and three-loop result. This justifies the choice \i = M t 
as default value. 
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